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ABSTRACT 

We develop the halo model of large-scale structure to include triaxial dark matter 
haloes and their intrinsic alignments. As a direct application we derive general expres- 
sions for the two-point correlation function and the power spectrum. We then focus 
on the power spectrum and numerically solve the general expressions for two different 
models of the triaxial profiles. The first is a toy-model that allows us to isolate the 
dependence of clustering on halo shape alone and the second is the more realistic pro- 
file model of Jing & Suto (2002). In both cases, we find that the effect of triaxiality is 
manifest as a suppression of power at the level of ~ 5% on scales k ~ 1 — 10 ft.Mpc _1 , 
which in real space corresponds to the virial radii of clusters. When considered by 
mass, we find that for the first model the effects are again apparent as a suppresion 
of power and that they are more significant for the high mass haloes. For the Jing 
& Suto model, we find a suppression of power on large scales followed by a sharp 
amplification on small scales at the level of ~ 10 — 15%. Interstingly, when averaged 
over the entire mass function this amplification effect is surpressed. We also find for 
the 1-Halo term on scales k < 10/iMpc -1 , that the power is dominated by ellipsoidal 
haloes with semi-minor to semi-major axis ratios a/c < 0.7. 

One of the important features of our formalism is that it allows for the self- 
consistent inclusion of the intrinsic alignments of haloes. The alignments are specified 
through the correlation function of halo seeds. We develop a useful toy model for 
this and then make estimates of the alignment contribution to the power spectrum. 
Further, through consideration of the (artificial) case where all haloes are perfectly 
aligned, we calculate the maximum possible contribution to the clustering. We find the 
hard limit of < 10%. Subject to further scrutiny, the proposed toy-model may serve 
as a means for linking the actual observed intrinsic alignments of galaxies to physical 
quantities of interest. 

Key words: Cosmology: theory - large scale structure of Universe - Galaxies: grav- 
itational clustering 



1 INTRODUCTION 

It is well known that the correlation functions for the dark 
matter contain a wealth of information about the cosmolog- 
ical model, the physics of the dark matter and the relative 
mix of baryons, photons and neutrinos (see Bond & Efsta- 
tiou 1984; Seljak & Zaldarriaga 1996; Eisenstein & Hu 1998). 
However, it is less well known that they also contain infor- 
mation concerning the shapes of the structures that form. 
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Under the assumption that the initial density field is Gaus- 
sian random, then during the linear stages of gravitational 
collapse all of the correlation functions, except for the 2-pt 
function, are zero and the fluctuations are effectively shape- 
less. However, under continued gravitational collapse, asym- 
metries develop and are amplified. For the hierarchical Cold 
Dark Matter (CDM) models this results in halo formation 
and then the emergence of the 'cosmic web' (Bond & Myers 
1996). Subsequently, all of the higher order correlation func- 
tions are non-zero and must now contain information about 
the morphological structure of the density field. 

Recent attempts to model the 2-pt correlation func- 
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tion have shown that qualitatively this can be matched 
with no further assumption than: that the haloes that form 
are spherical with some particular density profile (Navarro, 
Frenk & White 1997, hereafter NFW); that the large-scale 
over-densities are ellipsoidal (Sheth, Mo & Tormen 2000); 
and that halo formation is preferential in the regions of 
large-scale over-density (Sheth & Tormen 1999, hereafter 
ST). However, when this approach is used to model the 3- 
pt correlation function, significant discrepancies are found 
(Scoccimarro et al. 2001). Currently, the root of these dis- 
crepancies is believed to be simply due to the triaxiality of 
dark matter haloes. However, this has yet to be explicitly 
shown. 

Furthermore, whilst the 2-pt statistics can be qualita- 
tively reproduced with simplified models, it still remains to 
be seen whether these approaches can be made to precisely 
match results in general (Smith et al. 2003; Huffenberger & 
Seljak 2003). With current galaxy redshift surveys, such as 
the 2-degree Field Redshift Survey (Colless et al. 2001) and 
the Sloan Digital Sky Survey (Strauss et al. 2002) , being suf- 
ficiently large enough to produce hi-fidelity measurements of 
the 2-pt (Percival et al. 2001; Hawkins et al. 2003; Zehavi 
et al. 2002; Tegmark et al. 2004) and 3-pt galaxy clustering 
statistics (Jing & Borner 2004; Kayo et al. 2004), and with 
planned cosmic weak shear surveys expected to be capable 
of measuring the projected 2-pt and 3-pt matter clustering 
statistics to similar levels of accuracy (Aldering et al. 2004), 
resolving these issues becomes important. 

In this and subsequent work, we will explore these ques- 
tions in detail. Our first aim, in this paper, is to develop a 
self-consistent analytic model that allows us to make predic- 
tions for the nonlinear dark matter clustering signal, given 
information about the cosmological model and now, addi- 
tionally, the shapes of dark matter haloes and also their 
intrinsic alignments. A second aim will be to provide, as 
an application of the model, predictions for the 2-pt clus- 
tering statistics. We achieve these goals by developing the 
halo-model of large-scale structure (Seljak 2000; Peacock & 
Smith 2000; Ma & Fry 2000; and for a review see Cooray & 
Sheth 2003). 

An important by-product of this work is that the inclu- 
sion of halo shapes allows one to study the intrinsic align- 
ments of dark matter haloes alone. This is achieved through 
the inclusion of the correlation function of halo axis direction 
vectors, and we develop a toy-model to explore the effects 
of halo alignment on the power spectrum. 

The outline of this paper is as follows. In Section 2.1, 
we define some useful theoretical notions. In Section 2.2 
we present the triaxial halo model formalism and provide a 
derivation for the 2-pt correlation function. In Sections 2.3- 
2.5 we derive the power spectrum. In Section 3, we flesh-out 
the necessary details for performing calculations. In particu- 
lar we consider two interesting models for the density profile 
of triaxial haloes: the first is a toy model that we have devel- 
oped to explore how halo shape alone affects the clustering 
statistics; the second is the more realistic model of Jing & 
Suto (2002, hereafter JS02). Here, we also develop our toy 
model for the intrinsic alignments correlation function. In 
Section 4, we present our results, and finally in Section 5 we 
discuss our findings and present the conclusions. 

Throughout, we have assumed that the cosmological 
model is the concordance model (Wang et al. 2000) and 



that the linear power spectrum is given by Efstathiou, Bond 
& White (1992) with normalization <j$ = 0.9 and shape pa- 
rameter r = 0.21. 



2 THEORY 

2.1 Basic definitions 

In what follows, we will seek to compute the lowest order 
clustering statistic of interest, that is the power spectrum of 
mass fluctuations, -P(k). This is defined by 

{S(k)S(k')) = (2n) 3 P(k)S D (k-k') , (1) 

where S(k) is the Fourier transform of the density fluctua- 
tion field <5(x) = p(x)/p — 1 and p is the background den- 
sity; angle brackets denote the ensemble average and 8 D is 
the Dirac delta function. The Power Spectrum is itself the 
Fourier transform of the real space 2-pt correlation function 
£(r), which can be similarly defined 



C(r) = <5(x)5(x + r)) 



(2) 



In the above definitions we have considered general 
anisotropic density fields. However, in cosmological appli- 
cations it is usual to assume statistical isotropy and homo- 
geneity of the Universe, hence these quantities become func- 
tions of scalar arguments only. In the work that follows, we 
will not make the above assumptions, since we are dealing 
with anisotropic dark matter haloes, but will show that for 
the power spectrum the scalar arguments that we desire will 
naturally emerge from our formalism. 

It will also prove useful to define some basic relations 
for the triaxial haloes (see Chandrasekhar 1969, for a full 
treatise). We start by considering an heterogeneous triaxial 
ellipsoid that has semi- axis lengths a, b and c, where a < b < 
c, and orthogonal principle axis vectors e a , et and e c . We 
define the triaxial coordinate system (it!, O, $) with respect 
to the principle axes of the halo, where e c is taken to be 
in the z-direction. The radial parameter R traces out thin 
iso-density shells, or homoeoids, and the parameters O and 
$ are the polar and azimuthal angles respectively. In this 
system the Cartesian components are 

x — -it! cos $ sin ; y = -R sin $ sin 9 ; z = it! cos O. (3) 

c c 

It is to be noted that the ellipsoidal angles differ from those 
of the spherical coordinate system. The parameter R can be 
related to the Cartesian coordinates and axis ratios through 

c 2 ~ a 2 + b 2 + c 2 ' ( ] 

The benefits of this choice of coordinate system are now ap- 
parent: if all the homoeoidal shells are concentric, and if we 
pick coordinates with the same axis ratios as the triaxial el- 
lipsoid, then the density run of the ellipsoid can be described 
by a single parameter: 



P (v) - p(R) 



(5) 



Furthermore, the mass enclosed within some iso-density cut- 
off scale R cut , can be obtained most simply by 



M 



J 



, 3 . . , ab 
a rp(r) — 47r— 



-Rcut 



dRR?p(R) 



(6) 
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where the ellipsoidal coordinates have allowed us to circum- 
vent the problem of evaluating complicated halo boundaries. 

It is now convenient to define what we mean by a halo: 
any object that has a volume averaged over-density 200 
times the background density is considered to be a grav- 
itationally bound halo of dark matter. This leads directly 
to the following relation between the mass, radius and axis 
ratios: 

M 200 = ^7ri& (J)200p. (7) 

The above definition was adopted in order to be consistent 
with the mass-function of ST, which we utilize in what fol- 
lows. 



2.2 The triaxial halo model 

In earlier implementations of the halo-model it was assumed 
that all of the matter in the Universe was contained within 
spherical dark matter haloes with some particular density 
profile and some distribution of mass (Seljak 2000; Peacock 
& Smith 2000; Ma & Fry 2000). We now re-develop the halo 
model formalism making the important modification that 
haloes are not in general spherical, but instead are more 
closely described by the family of triaxial ellipsoids. 

To start, we characterize each halo in terms of a set of 
stochastic variables that describe both shape and orienta- 
tion: a given halo of mass M will therefore have principle 
axis vectors e Q , e b and e c with axis ratios a/c and b/c. The 
density run of a particular halo is described by the func- 
tion p(r, M, £ , a) = M U(r, M, £ , a), where U is the mass 
normalized profile, and we have introduced the shorthand 
notation £ = (e a ,e(,,e c ) and a = (a,b,c). The density at 
any point can now be expressed simply as a sum over the N 
haloes that form the field 



p(r) = ^2M i U(r-x i ,M i ,£i,a i ), 



(8) 



where Xi denotes the position vector of the centre of mass of 
the i th halo. Following Scherrer & Bertschinger (1991), we 
can re-write the sum in equation (8) using the substitution 



-> J dyidMd£daY^S D (x-x t )S D (M - M t ) 

i ' i 

x8 D (£-£i)8 D , (9) 

which allows us to consider integrals over continuous vari- 
ables rather than a sum over discrete quantities. The integral 
over £ in the above expression indicates an integral over all 
possible orientations of the halo. The orientation of the halo 
frame can be specified relative to a fixed Cartesian basis set 
through the Euler angles. These represent successive rota- 
tions of the halo frame about the z-axis by a, the j/'-axis 
by /3, and the z"-axis by 7 (see Mathews & Walker 1970). 
Hence the process of averaging over all possible halo orienta- 
tions can be performed by integrating over all possible Euler 
angles. 

Lastly, to compute the ensemble averages we integrate 
over the joint probability density function for the N haloes 
that form the density field, and sum over the probabilities 
for obtaining the N haloes (see McClelland & Silk 1977, for 
a similar approach): 



r Ns 

(•••) = ^2p(N d \V) / J[dMi d 3 Xi da, d£i 

xp(Mi, .., M Nj , xi, .., xjv 3 , ai, .., a N] , £ u .., £ Nj \Nj) .(10) 

Provided the volume of space considered is large, then 
p(N\V) is very sharply spiked around N = nV ^> 1, where n 
is the mean number density of haloes. We restrict our study 
to this case only. Also, the ensemble average of equation (8) 
gives the mean density of the universe, p. 

To compute the 2-point correlation function of the tri- 
axial halo field we substitute equations (8) and (9) into (2) 
and compute the ensemble average according to equation 
(10). The result is a sum of two terms, the first accounts 
for the correlation between points in a single halo and the 
second accounts for the inter-halo correlation (hereafter the 
1-Halo and 2-Halo terms): 



£(r)=e 1H (r)+e 2H (r) 



(11) 



S ( r ) 



= |j J dMdxd£daLM 2 p(x,M,£,a) 

x U(n - x, M, £, a) U(r 2 - x, M, £, a) - 1 ; (12) 



AT 

£ 2H (r) = — / dMidM 2 dxidx2dfid£2daida 2 Mi M 2 



x U(n -xi,Mi,£i,ai){7(r 2 - x 2 , M 2 , £2, a 2 ) 
x p(xi,x 2 ,Mi,M 2 ,£i,£ 2 ,ai,a 2 ) - 1 , (13) 

where the separation vector r = ri — r 2 . 

The joint probability density function in the 2-Halo 
term may be written: 



p(l,2)=p(l)K2)[l + f(l,2)] , 



(14) 



where we have adopted the short-hand notation = 
p(xi, Mi, £1, ai), and where £ s is the seed-correlation func- 
tion, which describes the relationship between a particular 
halo's characteristics and those of all the other haloes. In 
the spherical halo model this would just be the halo-bias of 
Mo & White (1996), however in the triaxial model the func- 
tion is more complicated, with halo orientation vectors and 
shapes being influenced by those of neighbouring objects. 
We explore this in greater detail in Section 2.5. 

Neglecting for the time being £ s (l, 2), all that remains 
to arrive at the correlation function in the triaxial halo 
model is to deal with the joint probability density function 
for a single halo's characteristics: p(l). We assume that a 
halo's orientation, position and mass are independent ran- 
dom variables, and that the halo axis ratios are dependent 
on mass only. Hence we have: 

p(x,M,£,a) = i^p(£)p(a|M), (15) 

where n(M) is the halo mass function. In order to provide a 
uniform probability for the halo orientation on the sphere, 
the density function p{£ ) takes the form 

p(£)d£ = p(a, f3,j)dctdf3 d"/ , 

= ^ 2~ da rf ( cos ^ d ~f' ( 16 ) 
where the variables are restricted to the ranges: > a > 2n, 
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> (3 > it and > 7 > 2ir. On substituting equations (14) 
and (15) into (12) and (13) we find: 



£ 1H (r) = j dM dxd£ daM 2 n(M) p(£) p(a\M) 
x C7(n - x)t/(r 2 - x) ; 



(17) 



C 2H (r) 



dMi dM 2 dxi dx 2 d£\ dS 2 da! da 2 Mi M 2 



x n(M 1 )n{M 2 )p{£i)p(£ 2 )p{ai\M 1 )p(a 2 \M 2 ) 
xf (l,2)[/i(n-xi)t/ 2 (r 2 -x 2 ) ; (18) 

where we have suppressed the explicit functional dependence 
of the density profile U on all variables except the position 
vector and the subscripts refer to either halo one or two. 

2.3 The dark matter power spectrum 

The power spectrum of the triaxial halo field can be evalu- 
ated by taking the Fourier transform of expressions (17) and 
(18) to yield: 

P(k) = P 1H (k)+P 2H (k) ; (19) 



P 1H (k) = — dMd£daM 2 n(M)p(£)p(a\M) 



xU(k)U*(k) ■ 



P 2H (k) = — / dMi dM 2 d£x d£ 2 dai da 2 Mi M 2 



(20) 



x n(Mi)n(M 2 )p(f i)p(£ 2 )p(ai|Mi)p(a 2 |M 2 ) 
xP s (l,2)(7(k)C/*(k), (21) 

where U is the Fourier transform of the halo profile and 
P s (l, 2) is the Fourier transform of £ s (l, 2). 

2.4 The 1-Halo term 

In this Section we focus on the 1-Halo term of equation (20). 
In order to solve this equation in its present form we are re- 
quired to perform a 12-D numerical integration. Evaluation 
of the complete expression would therefore be both numeri- 
cally noisy and slow. However, we can use properties of the 
elliptical coordinate system defined in Section 2.1 to simplify 
the general result significantly. 

We begin by re- writing the 1-Halo term of equation (20) 
in the form 



(22) 



(23) 



P m (k) = J dMdaM 2 n{M)p{a\M)W{'k) , 
where the window function is defined to be 



W(k) 



d£p(£) 



dr U(r, £) expik • r 



We now observe the following symmetry: averaging over all 
halo orientations for a fixed vector k is equivalent to av- 
eraging over the direction vector k at fixed halo orienta- 
tion. Hence, we are able to replace the integral J d£p(£) — > 
J dk/47r. This means that W and hence P 1H are now simply 
functions of the scalar length k. Consider next the r inte- 
gral; if we transform to the ellipsoidal coordinates and recall 



from earlier that under such a transformation the halo pro- 
file is simply a function of iso-radius R and also that the 
iso-density surface at which the halo is truncated is simply 
some particular value of R = P cu t , then the window function 
becomes 



W{k) 



dk 

4% 



dRR 2 U(R)^ I dR exp [ik • r(R)] 



.(24) 



Considering now the R-integral, if we write the Cartesian 
components of the r and k vectors in terms of ellipsoidal 
and spherical coordinates as 

r = R ^ cos $ sin 6, ^ sin $ sin 9, cos 9^) , (25) 

k = k(cos(f>k sin 9k, sin 4>k sin 9k, cos 9k) , (26) 

then the integrals may be evaluated using the standard re- 
sults (Gradshteyn & Ryzhik 1994): 



d<f> exp(iw cos 4> + iv sin < 



27rJ() Ju 1 + v 2 



da; cos(u:r)Jo (v\/ 1 — x 2 ^ =j (^\/ u 2 + v 2 ^ , 



(27) 



(28) 



where Jo is a Bessel function of the first kind and jo is a 
spherical Bessel function. After a little algebra, we arrive 
at the final result for the 1-Halo component of the power 
spectrum: 



P 1H (k) 



where 



dM dan(M) M 2 p(a\M)W(k) , 



and where 



JdRR 2 ^ 



U(R) 4nio[kRf(9 k ,< 



(29) 



(30) 



■ cos 2 <f)k + -~ sin 2 4> k -(31) 



f 2 (9 k , <t>k) = cos 2 6> fc +sin 2 9 k 



It is informative to consider two limits of the above formula. 
Firstly, if haloes are spherical, o/c = b/c = 1, then we re- 
cover the standard expression for the 1-Halo term (Peacock 
& Smith 2000). Secondly, if we consider the limit kR <C 1, 
then the jo-term is approximately unity and we have, as re- 
quired, the power simply being related to the effective num- 
ber density of objects (Peacock & Smith 2000). 

Through these efforts we have thus reduced the dimen- 
sionality of our integral from 12-D to 6-D, and in Section 4 
we recommend two numerical methods for solving integrals 
of this type. 



2.5 The 2-Halo term 

We now turn our attention to the 2-Halo term. Again, our 
aim is to manipulate the general result of equation (21) into 
a form that is more amenable to direct computation. How- 
ever, we now have to deal with the added complication of 
possible intrinsic correlations in halo properties and this will 
require us to provide a detailed form for £ s (l, 2). 
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2.5.1 No alignments 

To begin with it will prove useful to consider the simplest 
possible case, namely that for which there are no correla- 
tions between the halo seeds other than those between posi- 
tion and mass. In this instance the seed correlation function 
reduces to that of the spherical halo model and we may write 
(Mo & White 1996, ST), 

P s (k,Mi,M 2 ,£i,£ 2 ,ai,a 2 ) = 6(Mi) b(M 2 )P L (fc) , (32) 

where P L is the linear dark matter power spectrum, and 
b(M) is the bias function for haloes of mass M. On inserting 
this into equation (21), and once again using the symmetry 
between averages over £ and k, we find that the 2-Halo term 
with no intrinsic alignments can be written 

1 2 



P^(fc) 



where 



dM Mn(M) b(M)U(k,M) 



P L (fc) (33) 



U{k,M) = — 



j dap(a\M) J dkU(k, M, a, £) . (34) 

Equation (33) is identical to the standard spherical halo 
model result (Seljak 2000). Thus we may understand the 2- 
Halo term, for this case, to be simply the correlation of the 
spherically averaged triaxial haloes. However, owing to the 
complication of determining the spherically averaged profiles 
and the associated virial radius, we provide an alternate ex- 
pression for the 2-Halo term that yields more easily to direct 
computation. Considering again equation (21), if we take the 
inverse Fourier transforms of the density profiles, transform 
from the spherical polar coordinates to the ellipsoidal, and 
interchange the integral over the halo Euler angles to an 
average over k, then we find 



P™{k) = P L (fc) 



-_ J dMMn(M) b(M) J dap(a\M) 



JdRR 2 U(R) J dkjo [kRf{6 k ,4>, 



01 



(35) 



2.5.2 With intrinsic alignments 

Although the alignment of galactic/halo spins is an old sub- 
ject, dating back to work by Hoyle (1949), it has received 
relatively little attention over the past fifty years. However, 
in recent times it has been the focus of much work. This 
renewed interest has mainly been driven by a need to under- 
stand the possible contamination that galaxy shape align- 
ments may contribute to cosmic weak shear surveys (Brown 
et al. 2002; Heymans et al. 2004). Despite this renewed in- 
terest, the detailed physics underlying the process still re- 
mains the subject of much debate (Lee & Pen 2001; Mackey, 
White & Kamionkowski 2002; Porciani, Dekel & Hoffman 
2002), and as such no universally accepted model has yet 
been identified (Heymans et al. 2004). Moreover, these pro- 
posed models are not specified in a form that allows for 
direct insertion into our formalism. Alternatively, a number 
of authors have explored the alignment problem through di- 
rect numerical simulation. Croft & Metzler (2000); Heavens, 
Refregier & Heymans (2000); Jing (2002) have explored the 
correlation of halo ellipticities projected onto the sky and 
Hatton & Ninin (2001); Hopkins, Bahcall & Bode (2004) 



have studied the correlation of the scalar product of the 
semi-major axis vectors with separation. However, no sim- 
ple phenomenological model for the alignment correlation 
function, in the form that we require, exists. Thus in order 
to proceed our strategy will be to develop a toy-model for 
the intrinsic alignment correlation function C(|e c • e c |) that 
will allow us to probe, in a qualitative sense, the effects of 
halo alignments on clustering. 

We now outline some of the assumptions upon which 
the model is based since these will be useful in what follows, 
but we reserve the complete details to Section 3. Firstly, we 
assume that the intrinsic alignment between two haloes is 
simply manifest as a correlation of the semi-major axis vec- 
tors e c ; the e b and e a axes we will assume are uniformly 
random. Secondly, we assume that the spatial correlation 
function of the halo centres is independent of the halo ori- 
entations. Thirdly, we assume that the axis ratios, a, are 
uncorrelated from one halo to another, and that these quan- 
tities depend only on the mass of each halo. Under these 
assumptions equation (14) reduces to 

p(l,2) = p(l)p(2){l+£(xi,x 2 |Mi,M 2 ) + 

C(m|Mi,M 2 )[1 + £(xi,x 2 |Mi,M 2 )]} , ( 36 ) 

where n = ej • e 2 . On inserting the above expression into 
equation (13) for the 2-Halo correlation function, we find 
that the first term in the curly brackets gives the mean 
density squared, the second gives the 2-Halo term without 
alignments, the third integrates to zero through the integral 
constraint, J C(fi) dfi = and the fourth term represents 
the intrinsic alignment contribution to the 2-Halo term. This 
last term can be written, 



Cia (r) 



= — ■ / dM 1 dM 2 d£id£ 2 da 1 da 2 dx 1 dx 2 Mi M 2 
P 2 J 

x n(Mi) n(M 2 ) p(ai\Mi) p(a 2 \M 2 ) p(£i) p(£ 2 ) 
xUi(n -xi)[/ 2 (r 2 -x 2 )f (r|Mi,M 2 ) 
CV|Mi,M 2 ). (37) 

Fourier transforming the above expression leads us to the 
intrinsic alignment contribution to the power spectrum, 



Pl 2 A H W 



P L (fc) 



P 2 



dM 1 dM 2 da 1 da 2 MiM 2 n(Mi) n(M 2 ) 
xb(M 1 )b(M 2 )p(a 1 \M)p(a 2 \M)F(k) , (38) 



where, 



P(k) 



f d£idt 



d£ r f/i(k,£i)^ 2 (k,£ 2 )C(M|Mi,M 2 ) , (39) 



and where we have made use of the halo bias relation given 
by equation (32). Focusing on equation (39), again we can 
swap the integral over £ i to an average over k. However, 
owing to the alignment correlation function we must now 
specify the fixed coordinate system of £i, and we let this be 
the usual Cartesian system. Next, if we expand £ 2 in terms 
of the three Euler angles a, f3 and 7, where these are now 
rotation angles relative to the £ 1 system, and on realizing 
that el ■ e 2 — n = cos /3, then we have that 



F{k): 



dk_ dad^dn ^ (k) - a {Kaj/3 ^ )CMMuM2)j (40) 
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where the function 

tr 2 (k,a,/?,7)= ^ J d 3 i?'!7 2 (i?')exp[ik'(Q,/3,7)-r'(R')] 

= ^ JdR'R' 2 U2(R')4ir} [kR'f(e' ky cj>' k )] y (41) 

where r' is a vector in the £2 frame and where k' is the 
fc-vector rotated from the £ 1 frame into the £2 frame: 

k' = TZ(a, (3, 7 )k = TZ Z „ (j)TZ y , (/3)TZ Z (a)k . (42) 

1Z is the z — y' — z" rotation matrix (see Appendix A for 
the explicit form) and 9' k and <f>' k are the spherical polar an- 
gles of the fc'-vector in the £2 frame. Finally, on inserting 
equation (41) into (40) and taking the inverse Fourier trans- 
form of J7i(k), we arrive at the final answer for the intrinsic 
alignments contribution to the power spectrum 

Pia = f dM 1 dM 2 M 1 M 2 n(M 1 )n{M2)b{M 1 )b{M2) 



2np 2 

x / da 1 da 2 P(a 1 |M 1 )P(a 2 |M 2 )^Ji^ 

J C l C 2 

x J dRdR' R 2 R' 2 U 1 (R)U2(R') 

x J dkio[kRf(e k ,4>k)]V(e k ,cj> k ) ; 
where we have defined the useful function 
V = J dadfid-yjo [kR'f{e' k ,4>'k)] C{ t i\M 1 ,M 2 ), 
and where for completeness, 



0' k — arccos 



s/3 (k x ca + k y sa) + c/3k z 



ks6' k 



(43) 



(44) 



(45) 



(46) 



rj = (c/3cacy — sas7)k x + (c/3sac7 + cas7)k y — s/3cjk z (47) 

where we adopt the short-hand notation s = sin and c = cos. 

The full power spectrum, incorporating all of the triax- 
ial halo effects is thus 



P(k) = P 1H (fe) + P^(k) + Pff(k) , 



(48) 



where the three terms on the right hand side are given by 
equations (29), (35) and (43) respectively. These equations 
together comprise the main results of this paper. In the next 
section, we summarize the specific model details that are 
required to make direct calculations of these quantities. 



3 CALCULATION DETAILS 

As mentioned earlier, we use the ST mass function and halo 
bias relations. Since these are now widely known, we have 
reserved all details to Appendix B. 



3.1 Models for triaxial density profiles 

We consider two models for the triaxial density profiles. The 
first is a toy model that we have developed in order to ex- 
plore how the clustering statistics are affected by halo shape 



alone. We refer to this as the 'Continuity' model. The second 
is the more realistic model of JS02. 

• The Continuity model: The main idea here is that 
triaxial density profiles can be generated from the spherical 
dark matter haloes with some specified density profile. To 
see this, consider distorting a spherical halo along two or- 
thogonal axes and let this deformation be done in such a way 
that the volume is preserved. If the spherical radial length r 
is then exchanged for the ellipsoidal radial length R, given 
in equation (4), where we use the axis ratios of the ellipsoid 
formed from the deformation of the sphere, then the density 
structure of the new triaxial halo is completely understood 
in terms of the parameters that specified the spherical halo 
and the axis ratios. In the following we will start from the 
spherical density profile model of NFW, 

Pcrit 8 C S (M) 



ps(r) 



y = r/r (M) 



(49) 



2/(1 + y) 2 ' 

where ro is the scale radius and <5J is the characteristic den- 
sity contrast (hereafter subscripts e and s represent the el- 
lipsoidal and spherical models, respectively). If we distort 
this following the recipe described above then we get the 
ellipsoidal halo profile 

Pent S c e {M) 



Pe(R) = 



y = R/Ro (M) , 



(50) 



y(l + ' 

where Rq{M) and 5^(M) are the ellipsoidal scale radius and 
characteristic density contrast. As is the case for the spher- 
ical profile, Ro(M) and 8l(M) are not independent, but are 
related through the continuity of mass. Hence, 

_ 200c|/3 

e log(l+Ce)-Ce/(l+C e ) ' K ' 

where c e = R200/R0 is the concentration parameter. One 
is then left with determining a single free parameter. We 
now realize that if the volume is preserved under the dis- 
tortion, then the physical density at the scale radius must 
be equivalent for the spherical and ellipsoidal haloes, hence 
p s {r ) = pe(Ro), which implies that 5%{M) = 8 C S (M). We 
therefore simply require a model for Sg to specify the triax- 
ial haloes. NFW found this was well fit by 



st = 3000(1 + z c y 



(52) 



where z c is the redshift of collapse, which can be determined 
from the arguments of Lacey & Cole (1993). Finally, the 
cut-off radius of the ellipsoidal haloes can be easily obtained 
from continuity of the mass of spherical and triaxial haloes 
and is given by 



li>t 



ab 



-1/3 



T200 ■ 



(53) 



This model would be complete if the masses of ST and 
NFW haloes were equivalent. However they are not. The 
NFW halo mass is defined to be where the volume averaged 
overdensity reaches 200p cr it, as opposed to 200p for the ST 
case. If one considers a halo that has a density run given 
by the NFW model, then the different mass definitions sim- 
ply correspond to two different truncation densitys. Thus 
one may determine the mapping between the two mass def- 
initions through using the NFW profile and again equating 
the physical density at the scale radius. This time the im- 
portance of using the physical density is that the density 
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contrast relation is now modified by a factor of O. Hence 
one finds 



<5g T — 5npw/^ • 

Using equation (51) we find the relation 



1 / cnfw ' 

Q \ CST 



\ 3 log(l + C S t) - C ST /(1 + C S T 

/ log(l + cnfw) — cnfw/ (1 + cnfw) 



(54) 



= 1 ,(55) 



where cnfw = 7*^ o W / r o and where cst = f20o/ r o- This can 
be solved numerically to give the mapping between the two 
different concentration parameters and thus the mapping 
between the masses. In practice, once we have established 
the mass transformation we fit a spline function to it and 
use the relation 



CST 



= /l M ST \ 

" v n Mnfw / 



1/3 



cnfw 



(56) 



to convert between the different concentration parameters. 

Lastly, in this model we do not specify the probability dis- 
tributions for the axis ratios. Instead, we simply take these 
from the prescription of JS02, which is described below. 

• The JS02 model: A more realistic model for triaxial 
dark matter haloes was presented by JS02, who fitted an 
ellipsoidal NFW density profile to haloes measured directly 
from high resolution numerical simulations. Details of their 
model can be found in the original paper by JS02, but are 
nicely summarized in Oguri, Lee & Suto (2003). We discuss 
the points that are relevant for our application and refer the 
reader to these primary sources for a more complete account. 

The density run in the JS02 model has the same gen- 
eral form as the ellipsoidal NFW profile given in equation 
(50). However, now the normalization parameters S^(M) and 
Ro (M) are independent quantities and are to be determined 
as follows. Firstly, the virial radius of a spherical halo with 
equivalent mass to the triaxial halo is found. From this a 
scale length R e is determined. This is defined to be a fixed 
fraction of the spherical halo's virial radius, 



Re 



0.45r v 



(57) 



JS02 state that this corresponds to a scale where the volume 
averaged over-density within R e is given by 



5A V 



■<">(?) 



(58) 



where A v i r (fi) is the over-density of virialization from the 
spherical collapse model. The characteristic density is then 
found through 



<5e 



Ae Ct/S 



log(l + Ce) - C e /(1 + C e 



c e = R e /Ro(M) (59) 



where c e is the ellipsoidal concentration parameter. JS02 
then provide a separate model for c e . The independence of 
Ro and 5^ can now be seen, since although both parameters 
depend on c e , Ro also depends on R e , and 5° also depends 
on A e . Thus owing to the independence of Ro and 5%, for a 
given triaxial halo with axis ratios a/c and b/c, we no longer 
have a simple method for determining the cut-off radius for 
the haloes from the initial mass, since the volume averaged 
over-density is not preserved. Instead, we must determine 
.Rvir by numerically inverting the relation 



(60) 



M = A^RlSl [log(l + y) - y/(l +y)}; y = ^ 
c /to 

Note that again the definition of halo mass adopted by 
JS02 is different to that which we have adopted, where the 
JS02 mass is defined to be 



Mis = -Tvrl iT pA vil {n) 



(61) 



However, the mapping between the masses can be found 
through solving the relation 



/cjs_\ 3 A™ 
V c ST / 200 



log(l + c S t) - c S t/ (1 + c S t) 



log(l + Cjs) - C JS /(1 + Cjs) 



(62) 



where cjs = Rlf r /Ro- The mapping between the two con- 
centration parameters is then simply given by 



cst 



A vir (n) M ST 
200 M JS 



1/3 



cjs • 



(63) 



JS02 also provided a convenient fitting formulae for the 
probability density functions for the axis ratios a/c and a/6 
and these we use throughout. As a final note, JS02 advocate 
that c e be drawn from a log-normal distribution with mean 
concentration c e . Since the effects of a stochastic concentra- 
tion parameter have already been investigated in the halo 
model by Cooray & Hu (2001), we do not consider this here, 
and simply take c e = c e . 

3.2 Correlation function for halo alignments 

Following the discussion of intrinsic halo alignments in Sec- 
tion 2.5.2, we now describe the toy model for C(/x). We 
seek a function, defined on the range fi — [— 1, +1], that is 
symmetric about fi = that obeys the integral constraint, 
J dfj,C(fi) = 0, and that gives rise to a uniform distribution 
of \i when alignments are weak. The simplest functional form 
that will allow us to model these effects is to let C(a*) be lin- 
ear in Hence, 



C(m|Mi,M 2 ) = k{M 1 ,M 2 ) 



M-- 2 



; -1 < A» < 1 (64) 



where k describes the slope of the relation and is in general 
a function of Mi and M2. 

We now argue for how the mass dependence of k may 
arise and suggest a possible form for its scaling. Within the 
framework of the extended Press-Schechter model (Bond et 
al. 1991; Lacey & Cole 1993), one expects that at any given 
redshift, the highest mass haloes are more aligned than the 
lower-mass haloes. This can be understood on the grounds 
that since lower mass haloes form at higher redshifts they 
have more time to virialize and hence, particle orbits will 
become isotropic. On the other hand, higher mass haloes 
form at lower redshifts and therefore their shapes should still 
maintain some memory of the tidal fields that they collapsed 
within and this leads to a spin/shape alignment in nearby 
haloes. 

As a consistency argument to this picture, we remark 
that in the work of JS02 it is found that low-mass haloes are 
more spherical than high-mass haloes and therefore, since 
no particular alignment can be assigned to these objects, we 
should expect the correlation function to be smaller than 
for the higher-mass objects. A further piece of evidence sup- 
porting this view, comes from the numerical simulation work 
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of Hatton & Ninin (2001), who observed a positive mass- 
dependent trend in the correlation of angular momentum 
vectors of dark matter haloes with local large-scale struc- 
ture. While this is not exactly what we are considering here, 
it is at least suggestive of a similar effect. 

Finally, we will make the further assumption that the 
dominant alignment effect occurs between haloes of the same 
mass, and we will therefore introduce a Dirac delta function 
8 D (Mi — M2) into equation (64). One may argue for this 
as follows: it is more likely that haloes that collapse out of 
the same tidal field and at the same time will have their 
spin/shape vectors aligned. However, the veracity of this as- 
sumption has yet to be established. Nevertheless, the benefit 
of this is that the mass dependence of k can be reduced to 
a simple two — parameter function of M, which means that 
C(n) is also specified by just two parameters. 

With these ideas in mind, we take k to be a power-law 
in v = 8 CI it/a(M), where <5 cr it is the linear theory collapse 
density from the spherical model and a 2 (M) is the linear 
theory variance of mass fluctuations. Hence, 

k(^) = /to z/ , (65) 

where e > for high mass haloes to be more aligned and 
where «o sets the normalization for the M» haloes. Also, as 
v — ► 0, k — ► 0, and we get, as required, that C(fi) vanishes 
and the distribution of n is uniform over [— 1,+1]. Given 
the lack of detailed information concerning halo alignment 
correlations, it is difficult to assign values for e and Ko with 
any certainty. In what follows, we therefore simply explore 
effects for the parameter choices: e = 1.0,0.5,0.0 and kq = 
2.0, 1.0, 0.5. 

Before continuing, we note that more recent work by 
Jing (2002) and Hopkins, Bahcall & Bode (2004) has shown 
that there is indeed a mass dependence in the alignment of 
haloes and that it is broadly in agreement with the picture 
that we have described above. Furthermore, the studies by 
Hatton & Ninin (2001) and Hopkins, Bahcall & Bode (2004) 
have shown that the alignment effect decreases with radius 
as a power-law. We would now like to emphasize that this 
is exactly what our alignment model predicts. This is to be 
understood by the fact that the alignment correlation func- 
tion at zero spatial lag, which is given by equation (64), is 
modulated by the cluster correlation function (see equation 
36). Thus according to our model, the general alignment cor- 
relation function should simply be a scaled multiple of the 
cluster correlation function. We intend to fully explore this 
elsewhere. 



4 RESULTS 

4.1 Window functions & power spectra 

In this Section, we present results for the numerical integra- 
tion of the 1-Halo term's window function, given by equa- 
tion (30), and then the full power spectra given by equations 
(29), (35) and (43). However, before proceeding, we feel that 
it is necessary to mention the numerical algorithms we use 
to solve the equations. 

Although our analytical efforts have reduced the dimen- 
sionality of the various power spectrum integrals dramati- 
cally, from 12- to 6-D for P 1H and P^ A and from 18- to 12-D 



for P lA , the numerical integration of these equations using 
serial quadratures is extremely inefficient and intensely de- 
manding on cpu time. We therefore recommend the use of an 
efficient multi-dimensional integrator, such as the Korobov- 
Conroy algorithm (Korobov 1963; Conroy 1967) or the Sag- 
Szekeres algorithm (Sag & Szekeres 1964), both of which al- 
low fast and accurate evaluations of these expressions (tak- 
ing ~10 minutes to evaluate equation (43) on a standard 
2.8Ghz 32-bit Intel processor with heavy optimization af- 
forded through the Intel compiler, and where the number of 
evaluation points in the Sag-Szekeres algorithm was set to 
100 million). Both of the above algorithms have implemen- 
tations vended through the Numerical Algorithms Group, 
listed as routines: dOlfdf and dOlgcf . 

Figure 1 shows how, for the two different density pro- 
file models considered, the window functions, specified by 
equation (30) , change as haloes become more prolate. From 
inspecting the top panels of these figures, it is clear that 
the overall effect is small. We therefore take the ratio of 
the window functions with those for the equivalent spheri- 
cal model to inspect the effect more closely, where by equiv- 
alent spherical model we mean the window function for a 
particular triaxial profile model with a/c = b/c = 1. For the 
continuity profile model (bottom left panel), we find that as 
the prolaticity increases, the ratio gets smaller. This relative 
suppression in the window function is maximal on scales of 
the order the virial radius of the halo, being at most ~15% 
for the extreme, a/c = 0.4, prolate objects. Considering the 
haloes with JS02 profiles (right panels) , we find that the ef- 
fect is first seen as a small suppression on large scales and 
then as an amplification on small scales. We ascribe this am- 
plification to the fact that in the JS02 model 8^ increases as 
haloes become more ellipsoidal. 

Figure 2 shows the total power spectrum in the triax- 
ial halo model with no intrinsic alignments, obtained by the 
sum of equations (29) and (35). Note that here we plot the 
dimensionless power spectrum: A 2 (fc) = 47rfc 3 P(fc)/(27r) 3 
(Peacock 1999). Considering the triaxial haloes with the con- 
tinuity profile model (left panels), again the effect is weak. 
Looking at the ratio with the spherical case, we show that 
it is maximal on scales of the order the virial radius of clus- 
ters and is manifest as a suppression of power at the ~ 5% 
level. The oscillatory features that were seen in the window 
function have been washed out by averaging over the halo 
mass and axis ratio distributions. Considering the resultant 
power spectrum from the JS02 model (right panels) , we find 
a similar overall effect as for the continuity model, with a 
suppression of the order ~ 5% on scales k ~ 2hMpc~ 1 . 
The oscillatory features seen in the window function again 
have been smoothed out and the small scale amplification 
is no longer apparent. For the case where haloes are not in- 
trinsically aligned, we find differences between the spherical 
2-Halo term and the triaxial of the order a few percent for 
both models. However, these effects occur on scales where 
the 2-Halo power is sub-dominant to the 1-Halo power and 
therefore can be neglected. 

Figure 3, shows the contributions to the 1-Halo term 
from different mass ranges. For the continuity profile model 
(left panel), we find for each mass bin considered, that the 
power is suppressed relative to the spherical case, and that 
the effect is strongest for the high mass haloes, being at the 
level ~ 6%. There is also an small amplification of power 
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Figure 1. Variation of the window functions with increasing prolaticity. Left panels show effects for the continuity density profile model 
and right show the density profile model of JS02. fn the top panels, the solid line represents spherical haloes, and the dash, dot-dash 
and dotted lines represent haloes with a/c = b/c = 0.8, 0.6, 0.4 respectively. The bottom panels show the ratio of the window functions 
with the equivalent spherical haloes from each model and line styles have been preserved. In all cases the halo mass was set to be 
M = 1.0 x lO 15 /!" 1 !//©. 




0.1 1 10 

k [h/Mpc] 

Figure 2. Dimensionless power spectrum in the triaxial halo model. Left panels show effects for the continuity profile model and right 
show the JS02 profile model. Top panels show: the total power spectrum (thick solid line); which is the sum of P 1H (dotted lines) and 
Pj&A (dash lines); the total power from the spherical halo model (thin solid lines). Bottom panels show the ratios: total triaxial halo 
power to total spherical halo power (solid lines); triaxial 1-Halo to spherical 1-Halo (dotted lines); triaxial 2- Halo to spherical 2-Halo 
(dash lines). 
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k [h/Mpc] k [h/Mpc] 

Figure 3. Contribution to the triaxial power spectrum from different mass ranges. Left panels show the continuity profile and the 
right show the JS02 profile model. Top panels, the dash-line and the dotted lines, going from left to right, show the contribution 
to the 1-Halo term for the mass ranges: M > 10 14 /l _1 itf s ; W 14 h- 1 M Q > M > lO 13 fc- 1 M ; 10 13 /i _:l Ar© > M > W 12 h- 1 M Q ; 
10 12 /i -1 Mq > M > lO 11 ^ - 1 Mq. The thin dotted lines are the power spectrum split by mass for the equivalent spherical haloes. The 
thick dotted line is the 2-Halo term. Bottom panels show ratios of triaxial halo power to spherical for the same mass ranges as above, 
where we have adopted different and obvious line styles. 




k [h/Mpc] k [h/Mpc] 

Figure 4. Similar to Fig. 3, but this time showing the break down of the power spectrum with axis ratio a/c. In both panels, we 
show, going up from the dash line to the triple-dot dash line, the contribution to the 1-Halo term coming from haloes with axis ratios: 
1.0 > a/c > 0.85; 0.85 > a/c > 0.7; 0.7 > a/c > 0.55; 0.55 > a/c > 0.3. Again the left panel is the continuity model and the right is the 
JS02 model. 
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k [h/Mpc] k [h/Mpc] 

Figure 5. Contribution to the power spectrum from intrinsic alignments of the triaxial haloes, using the linear alignment model of 
equation (64). Left panel shows the variation with alignment normalization ko, and the right panel shows the variation with the power- 
law index e. 



on small scales. For the JS02 profile model (right panel), we 
find that there is a characteristic de-amplification, followed 
by a strong amplification of power as one goes from large 
to small scales. Interestingly, the large-scale suppression of 
power is stronger for the lower mass haloes, being of the or- 
der ~ 8% for the M ~ 5.0 x 10 n /i _1 M Q objects on scales 
k ~ 10/iMpc -1 . However, the small scale amplification effect 
is strongest for the high mass haloes, boosting the power to 
around ~ 10% for the M ~ 10 15 h^ 1 Mq on similar scales. 
Thus one sees that the largest effects of triaxiality are most 
apparent in the highest mass haloes. This result follows in 
accordance with the hierarchical picture of structure forma- 
tion: high mass haloes form at late times in the Universe and 
thus particle orbits have insufficient time to circularize and 
hence structures are more likely to be triaxial than spher- 
ical. This effect is built into the JS02 probability density 
function through a mass dependent scaling of the stochastic 
variable a/c, which operates so that higher mass haloes are 
more likely to be triaxial than the lower mass haloes (see 
equations 16 and 17 in JS02). This indicates that if one is 
purly intersted in the clustering properties of extreme mass 
objects, then halo triaxiality plays a more significant role. 

Figure 4 shows the contribution to the 1-Halo term from 
different ranges of the axis ratio a/c. Somewhat interest- 
ingly, in both cases, for k < 10/iMpc _1 the power spectrum 
is entirely dominated by the most triaxial haloes, a/c < 0.7, 
and that the spherical haloes contribute very little to the 
overall power. This observed effect: highly triaxial haloes 
contributing the most on large scales and spherical haloes 
accounting for more on small scales, can be understood en- 
tirely from the mass dependent scaling of the distribuion 
functions for a/c. 

Figure 5 shows the contributions to the power spectrum 
from the 2-Halo intrinsic alignment term. We find that, for 
the linear correlation model specified by equation (64), the 
intrinsic alignments contribution to the power spectrum is, 
~ 10 orders of magnitude smaller than the 2-Halo term with 



no alignmnets. As will become apparent from the following 
sub-section, this result is actually insensitive to our choices 
for ko and e. Considering the variation of the clustering with 
Ko (left panel), we find that increasing/decreasing ko simply 
uniformally increases/decreases P{a- Next, considering the 
variation of the power with the power-law index e (right 
panel), we find that as e increases/decreases the alignment 
power spectrum decreases/increases, and that the variations 
are very similar to those with Ko . This may mean that e and 
ko are not independent parameters. 

4.2 Maximally alligned haloes 

As a corollory to this section we discuss the limiting case of 
a maximally aligned halo distribution, that is, one in which 
the orientation vector of every halo is perfectly alligned with 
the others. Clearly, this scenario is unphysical, however it 
does allow us to place firm constraints on the largest pos- 
sible contribution that alignments may provide. In order to 
achieve this we simply set the correlation function to be 

C(p) = (5 D (/i± 1) - 1 ; -l<£t<l. (66) 

In this limit we find that the function T> in the intrinsic 
alignment term given by (equation 43) becomes 

V{e k ,4> k ) = 4n J d^io[kR'f(e' k \<t>k)] 

-J dad(cos/3)d 7 jo [kR'f(6' k ,4>' k )] , (67) 

where k" = TZ(0, 0, tp)\t and k' = 1Z(a, /3, 7)k, and where 
we have used the fact that for /i = ±1, we may combine 
rotations over a and 7 into a new rotation over ip = a — 7. 

Figure 6 shows the maximum possible contribution to 
the power spectrum from the intrinsic alignments of dark 
matter haloes that have the continuity density profile. We 
see that for this case, the alignment contribution has signif- 
icantly increased by roughly ~ 8 — 10 orders of magnitude. 
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Figure 6. Maximal contribution to the power spectrum from 
intrinsic halo alignments. The top panel shows: the maximum 
alignments (cross hatched region defines the result ± the rms 
error); the 2-Halo term with no alignments (dash line); the 1-Halo 
term (dotted line); and the total contribution from all terms (thin 
solid line). The bottom panel shows the ratio of the maximum 
to the total power. The triaxial density profile model was the 
continuity model. 

However, relative to the total power (bottom panel), it is still 
a fairly small quantity, being at most an effect of the order 
~ 10% on scales k ~ 1/tMpc -1 . Given that this alignment 
correlation function is physically unrealistic, we may now 
conclude that the true contribution to the power spectrum 
from alignments must be < 10%. 



5 DISCUSSION & CONCLUSIONS 

In this paper we have developed a formalism for calculat- 
ing the clustering statistics of a distribution of intrinsically 
aligned, triaxial, dark matter haloes: the triaxial halo model. 
This formalism facilitates the exploration of the importance 
of halo alignments and shapes on the density clustering 
statistics. 

As a direct application of the formalism, we have con- 
sidered the lowest order clustering statistic of interest, that 
is the power spectrum. We found, as usual, that the general 
result separates into a term that describes the clustering 
between haloes (2-Halo) and a term that describes the clus- 
tering within a halo (1-Halo). However, with the inclusion 
of intrinsic halo alignments we found that the 2-Halo term 
itself separates into two terms: the first described the halo- 



halo clustering without alignments; the second described the 
clustering due to halo alignments. We have derived compact 
analytic forms for these general relations, in all cases dra- 
matically reducing the dimensionality of the integral equa- 
tions. These were then solved numerically using two different 
multi-dimensional integrators. We considered two different 
density profile models. The first allowed us to explore pure 
shape effects on the clustering, the second was the more re- 
alistic profile model of JS02, which modified the halo central 
density based on halo shape. In both cases, we found that 
the effects on the power spectrum were, as expected, small, 
being at most ~ 5% for k ~ 1 — 10ft Mpc~ . However, when 
considerd by mass we found a more significant effect for the 
high mass objects, with a suppression of the order ~ 6% for 
the continuity model and an ~ 5% suppression of power on 
large scales followed by a ~ 10— 15% amplification of power 
on small-scales for the JS02 model. The effects of halo tri- 
axiality will be important to account for when interpreting 
precision measurements of cluster correlation functions on 
small scales. We have also found that the 1-Halo power is 
dominated by haloes with a/c < 0.7. 

We have explored the impact of halo alignments on the 
power spectrum. In order to achieve this it was necessary 
to develop a toy-model for the correlation function of the 
semi-major axis direction vectors of haloes. This model was 
constructed so that high mass haloes were more likely to 
be alligned than lower mass haloes. We further made the 
assumption that only haloes of the same mass are alligned. 
Although questionable, this assumption provided a means to 
consider the effects of a mass dependent alignmnet correla- 
tion function. For this model, we found that halo alignment 
contribution to the power spectrum was, surprisingly, ~ 10 
orders of magnitude smaller than the 2-Halo term without 
alignments. Modifying the normalization parameters for the 
C(n) does not significantly increase the alignment power. 
Thus, if our alignment model is correct halo alignments are 
completely unimportant for density clustering statistics. 

We then constructed the maximal alignment correlation 
function. This allowed us to explore the physically unreal- 
istic case where all haloes of all masses are completely al- 
ligned with one another. We found, for the triaxial haloes 
with the continuity profile, that the maximum alignment 
contribution to the power spectrum was ~ 10% on scales of 
the order k ~ 1 ftMpc -1 and significantly less on all other 
scales. This lead us to the conclusion that the true alignment 
contribution to the power must be < 10%. 

We now add an important and necessary cautionary 
note. In order to make accurate predictions for the clus- 
tering in the triaxial halo model, we require an accurate 
model for the statistical properties of the triaxial dark mat- 
ter haloes. Clearly, the first density profile model that we 
considered, the continuity model, was constructed simply 
as a toy-model, and as such the predictions should not be 
expected to match reality. In the second model, the JS02 
model, the profiles were constructed to match results from 
numerical simulations. However, we have discovered some 
aspects of this model that should be clarified before it can 
be used with confidence to make accurate predictions for the 
clustering statistics. 

Firstly, consider the special case of a halo that is ac- 
tually spherical and apply the JS02 formalism to it, having 
obtained the two independent normalization parameters, we 
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ask the question: At what radius does the average overden- 
sity reach that at which we defined the halo? The answer 
is not the same as the virial radius that we defined from 
the mass. Secondly, the ellipsoidal concentration parame- 
ter c e = R e /Ro depends on the axis ratio a/c, and not 
the second ratio b/c. This leads to the following problem: 
When averaging over the axis ratio distributions, we found 
that whilst a haloes shape and characteristc density may 
change dramatically through changes in b/c, the concen- 
tration paramter remains fixed. This results in inconsistant 
density structures for the haloes, since we found that haloes 
of a given mass and a/c could be less and more dense than 
spherical haloes, depending on the value of b/c. In reality 
it is more likely that haloes that are triaxial are either less 
dense, or of equivalent density, or more dense than spheri- 
cal haloes, but not all. Whilst the model of JS02 is ground 
breaking in many ways, we feel in light of these problems 
that some aspects should be re-visited. 

For the 1-Halo term we have found that modelling the 
density structure of haloes with triaxial ellipsoids produces 
an effect of the order 5 — 20%, relative to equivalently de- 
fined spherical haloes. For observational programs that hope 
to measure the small scale power spectrum to these levels 
of accuracy or better, one must therefore account for halo 
triaxiality when interpreting data with halo models. 

We mentioned earlier that if one takes halo concentra- 
tion to be a stochastic variable in the spherical halo model, 
then this too increases the small scale power (Cooray & 
Hu 2001). The effect is thus degenerate with triaxiality for 
the JS02 haloes. If the main reason for the stochasticity of 
the halo concentration could be attributed to the incorrect 
assumption that haloes are spherical, then one might un- 
derstand the degeneracy. However, JS02 found that there 
is a comparable scatter in the concentraion parameter for 
the ellipsoidal haloes as there is for the spherical. This they 
suggest is due to differences in the merger histories of the 
haloes and not asphericity. 

A further effect on the small scale clustering that we 
have not yet mentioned is that of halo subtructures. The halo 
model formalism was extended to include this by Sheth & 
Jain (2003). Effects on the density power spectrum were then 
explored by Dolney, Jain & Takada (2004). They found that, 
in general, the small-scale clustering becomes dominated by 
substructures below a certain scale. The scale depends on 
the sub-structure mass function and density distribuiton. 
Again, the effect appears to be degenerate with the effects 
of triaxiality. 

In this paper we have considered only the lowest order 
clustering statistics, the 2-pt auto-correalation function and 
the power spectrum. However, the overall goal of this work 
is to explore how the shapes of dark matter haloes and their 
intrinsic alignments influence the hierarchy of correlation 
functions. Since the 2-pt clustering statistics are the least 
sensative to the shapes of structures, it is not surprising that 
the degeneracies, noted above, have been found. However, it 
is expected that these will be broken through consideration 
of the higher order statistics, such as the bispectrum, and 
we will focus on this in a subsequent paper. 
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and where v' = y/av, a = 0.707, q = 0.3 and A w 0.322. 
Sheth, Mo & Tormen (2000) argued that the success of the 
ST formula was due to the fact that the ellipsoidal model for 
halo collapse was more realistic, as opposed to the spherical. 

We will also require a model for the halo biasing. Here 
we again look to the work of ST, who found 



APPENDIX A: ROTATION MATRIX 



Owing to there being several equivalent ways to define the 
Euler angles for the rotation matrix H(a,(3,j), we make 
explicit the definition that we use throughout. The matrix 
for the z — y' — z" rotation is given by: 



I c/3cac7 c/3sqc7 — s/3c7 \ 

— SQS7 +CQS7 

— c/3cqs7 — c/3sas7 s/3s7 

— SQC7 +CQC7 

\ s/3ca s/3sa c/3 



,(A1) 



where we have adopted the short hand notation c = cos and 



APPENDIX B: THE MASS FUNCTION AND 
HALO BIASING 

We model the mass function of dark matter haloes using the 
model of ST, since the predictions are in excellent agreement 
with the halo abundances measured directly from numerical 
simulations (ST Jenkins et al. 2001; Reed et al. 2003). The 
ST mass function is 



dn(M) 
dlogM 

where 



din v 



d\nM 



v = S col /a(M) 



*>=-(-i)(£W4 



(Bl) 



(B2) 
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